Singlet-triplet filtering and entanglement in a quantum dot structure 
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We consider two interacting electrons in a semiconductor quantum dot structure wliich consists of a small 
dot within a larger dot, and demonstrate a singlet-triplet filtering mechanism which involves spin-dependent 
resonances and can generate entanglement. By studying the exact time evolution of singlet and triplet states we 
show how the degree of both filtering and spin entanglement can be tuned using a time-dependent gate voltage. 

PACS numbers: 



In semiconductor quantum dots which can be formed for 
example at the interface of a gated GaAs/AlGaAs heterostruc- 
ture, we can trap electrons and control their number pre- 
cisely down to one''^'^"*. The fact that the electrons in the 
dots are well-confined within a region of some nanometers 
and they are therefore isolated from the remaining environ- 
ment of the host material, offers not only the opportunity 
to explore atomic-like effects and electron correlation s^'^'^'^ 
in the solid state, but also to consider the quantum dots as 
promising candidates for nanoelectronic applications. For 
instance the proposal for realisation of spin and charge 
qubits in double quantum dots^ has attracted a lot of in- 
terest both theoretically^ -^i'"-"''^ and experimentally ' 3- 14, 15. 1 6 
and the demonstration of entanglement, which is a necessary 
ingredient for two-qubit gates, is one of the main goals. Of 
course many other systems which include quantum dots make 
use of the Loss and DiVincenzc^ exchange-energy mech- 
anism based on the Coulomb interaction to generate spin- 
entanglement, such as for example flying qubits in surface 
acoustic wave structures'^ and static-flying qubits in quantum 
wires'** '^. The singlet-triplet qubit in a double quantum dot 
is also a promising candidate for quantum computation. The 
two-electron singlet-triplet (Sz — 0) states constitute the two- 
level quantum system and as has been shown this scheme, 
which is also under experimental investigation^^, is efficient 
for universal quantum gates^"'^'. 

In this paper we propose and investigate a semiconductor 
dot structure which can induce a mechanism to filter singlet- 
triplet states and thus to generate maximal spin entanglement 
between two electrons. The operation is controlled by a time- 
dependent gate voltage which can tune the degree of filtering 
and entanglement. 

To be specific, we choose material parameters for GaAs 
and consider the quasi-one-dimensional quantum dot struc- 
ture which is shown in Fig. lHa), loaded with two electrons. 
It consists of a small open dot which is formed in the centre 
of a much larger dot. The small dot is such that it can bind 
only one electron and this means that the Coulomb interac- 
tion forces the second electron to occupy the region of the 
large dot. The latter has width L ~ 800 nm ~ SOa^, where 
a*g is the effective Bohr radius for GaAs. In this regime the 
two electrons are in the strong correlation regime^^'^^ i.e. the 
Coulomb interaction dominates over the kinetic energy and in 
the ground state the two electrons are localised in regions for 



which the electrostatic repulsion is minimised. Note that this 
is the case for both singlet and triplet states which give virtu- 
ally the same electron distribution and they are separated with 
a very small antifeiTomagnetic exchange energy J. 

To explicitly demonstrate and analyse these effects we have 
studied the two-electron problem with exact diagonalisation^i 
considering the parabolic -band Hamiltonian 



2m* dx' 



+ Vcixi,X2), (1) 



where m* — 0.067mo is the effective mass of the electrons 
for GaAs. We have modelled the quantum dot with the Gaus- 
sian confining potential Vd{x) ~ — Vo exp(— a:^/2/^), for - 
400 nm < x < 400 nm and Vd{x) — oo otherwise^^. The 
parameters Vo , lo determine the depth and width of the small 
dot respectively and they are chosen to give only a single 
bound energy level -1.5 meV. The Coulomb term is given 
by Vc{xi,X2) = q'^/Anereor, with r = ^ {x\ - X2Y + \l 
and tr = 13 the relative permittivity in GaAs. This simplified 
form of the Coulomb interaction assumes that all excitations 
take place in the x-direction, whereas in the y, z-directions 
the electrons occupy at all times the lowest transverse modes. 
To satisfy this assumption we choose for all the calculations 
Ac — 20 nm, which gives a physical confinement length in 
the transverse y, z directions much smaller than that in x. The 
two-electron time-independent Schrodinger equation is solved 
numerically by the configuration interaction method. The 
low-energy eigenstates consist of two closely spaced singlets 
and two triplets with a somewhat larger energy gap to higher- 
lying states. Fig. [Hb) shows the lowest triplet eigenstate 
$^(a;i , X2) (antisymmetric) and Fig. [Ttc) typical energy lev- 
els. The coiTesponding effective charge density of $^(2:1, X2) 
which is derived by integrating the two-electron charge den- 
sity over the spatial coordinates of one of the two electrons i.e. 
p^{x) = 2^ / |$o (2:, X )\'^dx is shown in Fig.fTfa). For this 
distribution one electron is bound in the small dot, whereas 
the second electron is at the left and right corners with equal 
probability due to the symmetric potential. The charge densi- 
ties associated with all states in the low-lying manifold are in 
fact virtually identical to that shown in Fig. (Ha). 

This can be understood from a simple Heitler-London pic- 
ture using one electron states derived from a Hartree ap- 
proximation in which we replace the trapped electron in 
the small dot with an effective Hartree potential Vh{x) — 
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FIG. 1 : (Color online) (a) Quantum dot confining potential and effec- 
tive charge density (in arbitrary units) of the lowest triplet eigenstate. 
(b) Contour plot of the lowest triplet eigenstate. (c) The two-lowest 
singlet and triplet pairs. 
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FIG. 2: (Color online) Effective charge density for (a) left and (b) 
right triplet states. The quantum dot confining potential is also 
shown. 



For the dynamics that we describe below it is important to 
introduce left and right states which are formed by combin- 
ing the two lowest eigenstates for singlet $f = + $f , 
$1 = and triplet = -<f>J + $f , <i>]J = ^'^ + ^J 

(unnormalised). For these states one electron is bound in the 
small dot whereas the second electron is localised to the cor- 
responding corner as we demonstrate in Fig.|2]for the triplet. 



JVc{x,x )\ipcix )\'^dx , where "tjjdx) describes the single 
electron state of the small dot. The second electron feels the 
approximate potential Vh{x) + Vd{x), which for the lowest 
two-electron states of interest acts as a double-well poten- 
tial. The two lowest single-electron states tp-{x), 
(bonding-antibonding) of this potential have a small energy 
splitting and peak in the left and right corners. Using the 
one electron Hartree states we may form the two lowest elec- 
tron states (independent of spin) ^o{xi,X2) = ip-{xi)ipc{x2) 
and ^i{xi,X2) — tp+{xi)-4!cix2) and similarly for the next 
higher pair Within the Hartree model this means that not 
only the lowest singlet, triplet eigenstates have a distribu- 
tion of the form of Fig. [TJa), but also the first excited eigen- 
states which indeed agrees with the results that we derive 
from the exact diagonalisation. Note that even though the 
Hartree model gives great insight into the two-electron ener- 
gies and distributions of the low-lying states, with remarkably 
little effort, the symmetric and antisymmetric combinations 
<^'^''^{xi,X2) = i$o(a;i,a;2) ± $o(a;2,a;i)]/%/2 give a poor 
approximation to the singlet-triplet splitting, with even the in- 
coiTect sign. This is little improved by using more accurate 
one-electron states, e.g. from Hartree-Fock or DFT, being a 
consequence of strong correlations requiring accurate solu- 
tions of the two-electron problem. However, as we discuss 
later, the system may be described accurately by an extended 
Hubbard model. 



Although the effective charge density of the low-lying sin- 
glets and triplets is virtually identical, the wave functions 
are quite different and give rise to effective interactions be- 
tween spins, responsible for their entanglement. To demon- 
strate how these states can generate spin-entanglement and 
give rise to singlet-triplet filtering, we first prepare the sys- 
tem such that one electron is located on the left with spin up 
and the other in the central dot with spin down. A suitable 
non-entangled state of this type may be written as a super- 
position of the left singlet and the left Sz ~ triplet state 
i.e. '^^i = ($f xfi + ^'[xji) / V^' with the spin eigenstates 

xff (1,2) = [xT(l)Xi(2) TXi(l)XT(2)]/V2. An approxi- 
mation to this state may be obtained in principle by applying a 
small source-drain bias, Vsd, across the quantum dot structure 
and allowing the system to relax to its ground state (singlet). 
The spins may then be initialised using magnetic fields and 
a microwave pulse, as suggested for scalable qubit array*^. 
If Vsd is removed, then oscillations (in a similar fashion to a 
double dot system) and entanglement will develop with a typ- 
ical time for the triplet component of the wave function on 
the order of nh/AE'^ ^ 1 //s, where AE^ ^ Ej" - Ej is 
the energy splitting of the two lowest triplet eigenstates when 
Vsd=0 (similarly for the singlet). To speed-up the process we 
can apply a time-dependent gate voltage that will tune in such 
a way the dot potential so that to increase this energy splitting. 
In this work we have modelled the gate voltage potential using 
the expression Vg(x, i) = —Vp{t)e-Kp{—x'^/2lg), (Ig ^ 130 



nm) which is driven by a triangular pulse of the form 

f T^t, < < < Tp/2 

Vp{t) - < ' (2) 

[-^t + 2Vi„ Tp/2<t<Tp, 

where Tp is the period of the pulse and Vb the maximum gate 
voltage. We have studied the case for which the gating rate 
a = \dVp/dt\ is such that to a good approximation only the 
two lowest singlet and two lowest triplet states are involved in 
the dynamics. During the first half of the cycle (0 < t < Tp/2) 
the effect of the voltage is to decrease the effective width of 
the dot thereby increasing the energy splitting and the inter- 
action, whereas during the second half of the cycle {Tp/2 < t 
< Tp) the process is reversed. 

The dynamics of the two electrons is governed by the time- 
dependent Schrodinger equation with the Hamiltonian (1) and 
for a total potential Vd{x) + Vg{x, t) which is the sum of the 
dot confining potential and the time-dependent gate poten- 
tial. For t > 0, the spin eigenstates are unchanged for sin- 
glet or triplet components, whereas the evolution of the corre- 
sponding orbital states is given directly by the solution of the 
time-dependent Schrodinger equation which is implemented 
numerically using the staggered-time algorithm proposed by 
Visscher— . The time evolution of the initial state, is then 
determined by adding the separately determined singlet and 
triplet components. In Fig. |3] we show two examples of the 
final electron distribution for singlet and triplet components. 
Specifically, in Fig. 13 a) the singlet corresponds to a good ap- 
proximation to the right state, whereas the triplet to the left 
state. This is the ideal filtering regime (together with the op- 
posite limit) in which at the final time the two components 
occupy different spatial regions. Note that in this regime the 
two electrons are fully entangled provided the measurement 
domain is restricted either to the left or right region, detecting 
the singlet or the Sz — Q triplet state respectively, which are 
fully entangled states. Performing an additional measurement 
but initialising the electron spins to a purely triplet state (ei- 
ther of S'z = ±1) reveals the occupation region of the triplet 
state hence offering a way to discriminate between singlet and 
triplet components. Note that in the most general case which 
is shown in Fig.|3jb) the final two-electron state is a superpo- 
sition of left and right states for both singlet and triplet com- 
ponents. Clearly, a measurement with restriction to either the 
left domain or the right domain would probe a partially entan- 
gled state. 

To understand the dynamics and the difference between sin- 
glet and triplet components we can write the time-dependence, 
for example of the singlet, using the corresponding two low- 
est energy states at all times i.e. <i>'^ = Cf + C|$|, in 

the left-right basis with the amplitudes Cf = e^'^f costi;'^ 
and — ie~^'^^ sinw^. The parameters (f^ ~ JjJ[_Ef (i ) + 

(t')]d//2Sandw^ = ll\Ef(t)~El(t)\dt/2hA&^&XiA 
on the two lowest energy levels and change with time. Note 
that Lo^ determines the amount of filtering because it deter- 
mines the period of the oscillations between left and right 
states, whereas both Lp^ and uj^ determine the degree of en- 
tanglement as we demonstrate below. In Fig. |4] we show the 
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FIG. 3: Effective charge density for singlet and triplet components at 
the final time for (a) a case in the filtering regime and (b) the most 
general case. 
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FIG. 4: Energy splitting of the two-lowest levels for singlet and 
triplet components, as a function of time for one cycle in the filtering 
regime. 



energy splitting of the two lowest levels (the energies are cal- 
culated from instantaneous solutions) as a function of time 
for one cycle in the filtering regime which gives fully entan- 
gled electrons in time Tp ^0.3 ns (for Vf, ^2.8 meV). We see 
that the splitting is larger for the singlet and this means that 
uj^ > lli^ i.e. the singlet component oscillates faster than the 
triplet. We can understand this effect with a Hartree approx- 
imation: because the small dot has only one bound energy 
level, £c, the tunneling time (the time that the electron in the 
left corner needs in order to move to the right corner) for the 
triplet is longer due to the Pauli blocking. Within the Hartree 
model the electron in the large dot needs to tunnel through an 
effective double barrier due to the Coulomb repulsion from 




FIG. 5: (Color online) Dependence (at the final time), as a function 
of maximum gate voltage for fixed gating rate, of: (a) absolute value 
of cos (5 and right probabilities for singlet and triplet components; 
(b) concurrence, spin-flip and non spin-flip probabilities calculated 
for the right region. 



FIG. 6: (Color online) Dependence (at the final time) of quantities in 
Fig.[5]as a function of pulse duration time for two different maximum 
gate voltages, (a) 14 =2.75 meV and (b) 14 =5 meV. 



the trapped electron in the small dot, which has only one res- 
onance energy level ~ + C/c (with Uc the Coulomb energy 
when both electrons occupy Ec) and corresponds to a singlet 
state. For a triplet resonance to exist the small dot needs to 
have at least two bound levels due to the Pauli principle. It 
is worth noting that the efficiency of this singlet-triplet filter- 
ing mechanism and the induced entanglement has been also 
demonstrated in scattering problems between a static and a 
flying qubit'**''''. 

To quantify the degree of entanglement which arises as a 
consequence of the interaction when the electrons are close 
together we have calculated concurrence^^ at the final time, 
from spin-flip, P,/, and non spin-flip, P,isf, probabilities. For 
example at the right region the amplitudes for these probabil- 
ities are 



Cn. 



of 



c. 



(3) 



which are derived by writing the corresponding final state 
(ty) in the basis of spin eigenstates. Concurrence at the 
right region is then given by 



Psf+Pnsf 

with Pnsf = \Cnsff and P^j = IC^/P, thus 

m + Pir-^PiPK cos' S]^^ 



pS , pT 



(4) 



(5) 



where we have set P^ — |C|! p = sin^ w"^ and P]^ = 
|C^P = sin^ cj-^ and the relative phase 5 = Lp^ — (p^ . A 



similar expression can be derived for the left and even for the 
total region. Note that by definition < Cr < 1 where the 
limit Cu ~ corresponds to unentangled electrons, whereas 
the limit C^t = 1 to fully entangled electrons. The measure- 
ment in the right region is meaningful only when P^ ^ 0, 
and/or P^ ^ i.e. when the right state is occupied for sin- 
glet and/or triplet. We see from Eq. (|4| that Cr — 1 when 
Psf=Pnsf or equivalently from Eq. ^ when cos S — and/or 
when P^ — and simultaneously P^ ^ or vice versa. 

In Fig.|5]we present the dependence of probabilities, phase 
shift and concurrence as a function of maximum gate voltage 
Vb for a fixed gating rate a 18 meV/ns which ensures that to 
a good approximation only the two lowest eigenstates for both 
singlet and triplet components are involved into the dynamics. 
Figure |5j a) shows the probabihties Pj, P^ and the absolute 
magnitude of the quantity cos 5 calculated at the final time, 
versus Vb- We see that the probabilities P^, P^ do not oscil- 
late for small gates voltages Vb ^2.5 meV (but rather they in- 
crease) because the energy splitting of the two-lowest levels in 
this limit remains relatively small and therefore the tunneling 
time long. Further increase of the gate voltage induces well- 
defined oscillations. In particular the probability of the singlet 
component oscillates faster than the triplet because the energy 
splitting which determines the frequency of the oscillations is 
larger for the singlet. Note also that the frequency of oscilla- 
tions increases with gate voltage for both components follow- 
ing the increase of the energy splitting. Figure|5lb) shows the 
concurrence Cr and the spin-flip Pgf and non spin-flip Pnsf 
probabilities calculated at the final time. When these proba- 
bilities are equal the entanglement is maximum whereas when 
one of them is zero the electrons remain unentangled. For a 
gate voltage Vb ^2.8 meV we have an example of the ideal fil- 
tering regime for which P^ ^0, P^ ^ 1 and Pgf = Pnsf giv- 
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FIG. 7: Illustration of the two-electron spin dynamics in a three-step 
process I, II, III. The process (a) is allowed only for singlets whereas 
the process (b) is allowed for both singlets and triplets. 



ing Cfl=l, since only a triplet state occupies the right region. 
The generic condition for ideal filtering is when sin u!^=l and 
simultaneously sin u!'^=Q or vice versa. Other filtering cases 
(non-ideal) which give Cr =1 occur when one of the compo- 
nents is zero, say =0 and the other is non-zero P]^ ^0, but 
not however equal to one. Finally an interesting limit is when 
P^ — P]^ for which the concurrence reduces to the simple 
expression Cr = \ smd\ which depends only on the relative 
phase between singlet and triplet components. 

Figure |6] presents the dependence of various quantities as 
a function of the pulse duration time for two different gate 
voltages. As we see by comparing Figs. |6|a) and |6lb) the 
period of the oscillations can be controlled with the value of 
the maximum gate voltage with the larger Vf, inducing faster 
oscillations due to the larger energy splitting. 

Finally, we point out that the low-lying spin states relevant 
to the dynamics of this entangler which involve virtually iden- 
tical charge distributions at the left, centre and right of the dot 
structure suggest that this charge-spin system may be mod- 
elled by a three-site Hubbard model of the form 

Heff =s{nL + Ur) + Ecnc 

cr,i=L,R (6) 

+ U{nL,^nL.,i + nR^^UR^i) + UcUc^Tnc^i 
+ V{nLnc + ncriR), 

where ni = J^a^ho- = Y.a'^i.a'^h'y' ^ = L,R,C and 



C] ^{Ci^cr) creates (destroys) an electron on site i with spin- 
cr. The on-site orbital energy for the left and right sites is e 
and for the central site is Sc- The on-site Coulomb energy is 
denoted by U for the left and right sites and by Uc for the 
central site. The nearest site Coulomb energy is V and finally 
7 expresses the hopping between nearest sites. Note that all 
the physical parameters (e, £c, U, Uc, V, 7) may be estimated 
from the Hartree approximation described earlier Working 
within the restricted subspace for which the two-electron ba- 
sis states consist of six singlets and three triplets, we can 
extract the correct energy splitting of the two-lowest eigen- 
states (for both singlet and triplet), the corresponding eigen- 
states and the antiferromagnetic exchange energy. Solution of 
this time-dependent model Eq. (|6]) does indeed show qualita- 
tively the same behaviour as the original continuous problem 
but with considerable saving in computer time once the time- 
dependent parameters are known. The analysis based on the 
effective Hamiltonian Eq. ^ indicates that the physical be- 
haviour that we have demonstrated can also be realised in a 
triple quantum dot structure. More importantly, the Hubbard 
model gives insight into the behaviour of the system since it 
is readily mapped onto an effective charge-spin model for the 
low-lying manifold of two singlets and two triplets. In par- 
ticular, for singlets there are two processes by which an elec- 
tron may tunnel from left to right, as shown in Fig. |2l The 
process (a), that is only allowed for singlets, has amplitude 
Ai = 7^/ {U — V + Ec — s) ^ J where J is the Heisenberg 
exchange energy between an electron spin on the central site 
and one on either the left or the right site. On the other hand, 
process (b) is valid for both singlets and triplets and has am- 
plitude A2 — -f^ / {e ~ Ec ~ V). The relative rate of tunneling 
for singlet and triplet is thus tuned by the gate by changing Ec, 
since making Ec more negative the energy denominator in Ai 
decreases whereas that of A2 increases. 

In summary we have presented a quantum dot structure and 
described the two-electron distribution for the lowest singlet 
and triplet states. By studying the electron dynamics due to a 
time-dependent gate voltage we showed how we can induce 
a singlet-triplet filtering based on the Pauli blocking effect 
for the triplet state. This can generate full spin-entanglement 
within the order of ~0.3 ns. Both the degree of filtering and 
entanglement can be efficiently controlled with the gating rate 
and the maximum applied gate voltage. 

G.G. thanks UK EPSRC for funding. This work was sup- 
ported by the UK Ministry of Defence. 
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